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Background 


A concept for "Gasdynamic Fusion Propulsion System for Space 
Exploration" (Reference 1) has been theoretically investigated at the 
University of Michigan for the Department of Energy. An experiment to 
determine the feasibility of this concept is currently being pursued at MSFC 
by Bill Emrich. The goal of this program is to develop a fusion propulsion 
system for a manned mission to the planet Mars. Several significant 
technological hurdles stand between this concept and flight hardware. 
Most, if not all, these hurdles for this concept stand on the problem of when 
the deuterium and tritium atoms under fusion, the resulting neutrons and 
alpha particles are emitted in random directions (isotropically). The 
probable direction of emission is equal for all directions. This results in 
most of the energy being wasted, massive shielding and cooling 
requirements, and serious problems with the physics of achieving fusion. 

These problems could be alleviated if the atoms could be induced to 
undergo fusion and emit the neutrons and alpha particles with a preferential 
direction. The neutrons and alpha particles are always emitted in exactly 
opposite directions due to conservation of momentum. If the nuclear 
magnetic spin moments of the deuterium and tritium nuclei could be 
precisely aligned at the moment of fusion, the stream of emitted neutrons 
could be directed out the rear of the spacecraft for thrust and the alpha 
particles directed forward into an electromagnet to produce electricity to 
continue operating the fusion engine. This would result in the shielding and 
radiant cooling requirements being reduced to practically zero. The only 
thrust exhaust from the spacecraft would be electrically neutral neutrons, 
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thus preventing charging of the spacecraft. Also, in this configuration, 
extraction of energy from the fusion process would be almost 100%. 

Alignment of the nuclear magnetic spins of atoms occurs every day in 
hospitals around the world. The very useful diagnostic technique, NMR 
(Nuclear Magnetic Resonance) uses a magnetic field to precisely align the 
orientation of atomic nuclei in a patients body. The interaction of the 
magnetic field of the NMR machine and the nuclear magnetic spin of the 
nuclei causes the affected nuclei to align their spins in a precise direction. 
Radio frequency energy is applied, causing the spins of the nuclei to reverse 
orientation, or "flip". A sensor detects this spin flip and precisely locates 
the position of the atom whose spin has flipped. 

The concept of using a magnetic field to align the nuclear magnetic 
spins of atoms is a common practice. A magnet with a magnetic field of the 
proper orientation could be located near the interaction region where the 
atoms undergo fusion in the spacecraft propulsion system. The nuclear 
magnetic spins of the deuterium and tritium nuclei would be oriented in the 
proper position to direct all the neutrons rearward as thrust and the alpha 
particles forward to generate electricity. Approximately 80% of the energy 
in fusion is released as the kinetic energy of the neutron, the other 20% of 
the energy is carried by the alpha particle. Less than 10% of this generated 
electrical power would be required to sustain continued operation of the 
fusion engine. If this technique could be applied to fusion propulsion, a 
manned mission to Mars could be a near-term possibility. 
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Introduction 


In the gasdynamic fusion propulsion concept described in Reference 
1, the 17.6 MeV of kinetic energy produced by the fusion of deuterium and 
tritium nuclei appears as a 14.1 MeV neutron and a 3.5 MeV alpha particle. 
Some of this energy manifests itself in the form of bremsstrahlung and 
synchrotron radiation, which eventually appears as heat and must be 
disposed of by a space radiator cooling system. The material mass of the 
required space radiator is estimated by dividing the total thermal power 
produced in the fusion engine by the power radiated per unit mass of the 
radiator. This is frequently assumed to be 5 MW/metric ton. 

For a fusion system described in Reference 1 which utilizes a 
superconducting magnet with a total current density of 250 MA/m , a 
deuterium and tritium plasma density of 10 17 cm" 3 at a temperature of 20 kev 
that exits through a 3 cm nozzle producing thrust, then the space radiator for 
such a spacecraft propulsion system would have a mass on the order of 
21,000 metric tons. The wasted excess neutron power in this example is 
about 1725 MW/m 2 and the surface heat flux due to the bremsstrahlung and 
synchrotron radiation power is about 30 MW/m 2 . 

A concept for such an interplanetary space vehicle appears in Figure 
1. This diagram clearly shows the massive size of the radiator required to 
eliminate excess heat produced when neutrons are emitted isotropically 
along the axis of the fusion propulsion engine. 


6 



_ 

Neutral beam 
injector 



Direct energy converter 



Fuel tanks Habitat 


Figure 1. Interplanetary Space Vehicle Concept (Reference 2) 
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Nuclear Magnetic Moments and Spin Precession in a Magnetic Field 


The nucleus of an atom has an angular momentum due to its spin. 
This angular momentum is quantified (Reference 3,4,5) by the nuclear spin 
quantum number I multiplied by Plank’s constant h divided by 27T, or Hf . 
In addition, the nucleus also possesses a magnetic moment. This nuclear 
magnetic moment is expressed in terms of a nuclear magneton M defined 
by the equation 

_ e*f\_ 

in which e is the charge of the nucleus and m n is the mass of the nucleus. 
The nuclear g factor relates the magnetic moment G of a nucleus to its spin 
angular momentum. The nuclear g factor is defined as the ratio of the 
nuclear magnetic moment, expressed in units of nuclear magnetons, to the 
spin angular momentum, expressed in units of h . Thus, 

- M 
3 ~ I 


And 

e ^ 

/VI — a I = $ ^ ^ *** 


When a nucleus of magnetic moment G is in a constant magnetic 
field of induction B , it will precess about the direction of B with a 
frequency given by Larmor’s theorem 
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/v\ B 

[y ~ 1 h 


The magnetic moment G of a nucleus is thus determined by the 
Larmor frequency \7 which the nucleus of spin quantum number I acquires 
in a known constant magnetic induction B . 

The measured nuclear magnetic dipole moment G for a deuterium 
nucleus is + 0.8574376 nuclear magnetons with a spin of + 1 and for a 
tritium nucleus it is + 2.978962 nuclear magnetons with a spin of + Vi. 

When an accelerated beam of deuterons and tritons enter a region 
with a uniform external magnetic field, the spins of the nuclei are aligned 
with the magnetic field and precess about the magnetic field lines at the 
Larmor frequency. The evolution of the spin direction (References 6,7,8) 
and the orbital motion of a polarized beam in external magnetic fields is 
governed by the Thomas-BMT equation 



and the Lorentz force equation: 






where the polarization vector P is expressed in the frame that moves with 
the particle. From comparing these two equations it can readily be seen 
that, in a pure vertical field, the spin rotates Gy times faster than the orbital 
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motion. Here G is the nuclear magnetic moment of the deuteron or triton 
and gamma is the charge to mass ratio. In this case the factor G v then gives 
the number of full spin precessions for every full revolution, a number 
which is also called the spin tune frequency. 

The acceleration of spin polarized beams of nuclei is complicated by 
the presence of numerous depolarizing resonances. During acceleration, a 
depolarizing resonance is crossed whenever the spin precession frequency 
equals the frequency with which spin-perturbing magnetic fields are 
encountered. There are two main types of depolarizing resonances 
corresponding to the possible sources of such fields: imperfection 
resonances, which are drive by magnet errors and misalignments, and 
intrinsic resonances, driven by the magnetic focusing fields. 

By introducing a 180 degree magnetic spin rotator of the spin about 
the horizontal axis, the stable spin direction remains unperturbed at all times 
as long as the spin rotation form the 180 degree spin rotator is much larger 
than the spin rotation due to the resonance driving fields. Therefore the 
beam polarization, is preserved during acceleration. Such a magnetic spin 
rotator can be constructed by using either solenoidal magnets for low 
energy (less than 1 MeV) beams or a sequence of interleaved horizontal and 
vertical dipole magnets producing only a local orbit distortion for high 
energy beams. 
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Nuclear Spin Quantum Mechanics 


In the absence of spin, the Hamiltonian for the motion of a deuteron 
or triton with a charge of +e in an electromagnetic field is (References 
9,10,11) 

H 0 = it (?'£*)' + 

where A is the vector potential and (f) the Coulomb potential. Since the 
spin interacts with the magnetic field, the nucleon gains additional potential 
energy. The nuclear magnetic moment reads 

/ - / el \ * 

/H = 3 (t^T- ) 5 -■ V* ' 

where /tg = (eh/^mc) and the potential energy in the magnetic field is: 

\J = -/V\ ■ D 

The Hamiltonian of a nucleon with spin takes the following form: 

H - H. + yA b 5 - • B 

Knowing the value of the Lande g factor, this Hamiltonian leads to the 
Schrodinger equation of a particle with spin, also known as the Pauli 
equation. This is the equation for a system of two coupled differential 
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equations for spin wave functions describing the nuclei with the z 
component of their spin either up or down. Because of the form of the Pauli 
spin matrices, the spin system described by the Schrodinger equation is 
decoupled for the z component of the spin and only coupled by the x and y 
components. This Pauli equation is given by: 






■b e <^> * yMt 



To calculate the current density which results from the Pauli equation, the 
equation is rewritten in the form: 





The adjoint to this equation reads: 

+ 


«.v w*-**) 


t 




(fz ■ B 


because & is Hermitian and the magnetic field B is real. Rearranging and 
manipulating the Pauli equation and its adjoint yields the continuity 
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in the form 


c?lV 

Jf 


- 


J 


^ II 

where w= T ^ is the probability density and the current density of the 


nuclei is given by 




j 


a*- 



y/ + v<? - (W f J vj- ^ /4 ^ V 


With both components of the spin wave function ( spin up and spin down) 
and their complex conjugates, the total current density equation is 



or, rearranging, 



W) - Si * 
K v ' f -)- Si* ^ <K 


The probability density and the current density are composed 
additively of the parts of the two different spin directions; multiplication of 
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the nucleon current density j by the nuclear charge yields the electrical 
current density j t . 

The current density j e does not contain the spin; rather it is the current 
density caused by the orbital motion of the nucleon, each with its own spin. 
However, the spin of a nucleon also causes a magnetic moments, which can 
be expressed by a corresponding current called the spin current density j s . 
This current density cannot occur in a continuity equation in which the 
charge conservation is expressed by convection currents. 

The spin current density can be calculated using Maxwell’s equation 
for the curl of the magnetic field: 

Curl 

In this case the magnetization is replaced by the averaged 

density of the magnetic moment AA , where averaging is performed over 
the spin states. The nuclear magnetic dipole density is given by: 

yCA a ty* a V 



The contribution to the spin current density causing the magnetic 
moment of the nucleons is then described by: 
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j, = -CyU t Curl Lf* + 2Y 
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Kinematics of Nuclear Reactions 


When a deuteron and triton collide with sufficient energy, a nuclear 
reaction occurs similar to inelastic scattering where the incident and target 
nuclei are transformed into different nuclei or particles ( a neutron and an 
alpha particle). Regardless of the circumstances of the reaction, four very 
important quantities (References 12 and 13) are observed to always be 
conservedj 1) mass-energy, 2) linear momentum, 3) angular momentum, 
and 4) spin. The last three are critical to the concept discussed in this 
report. 

Conservation of linear momentum is a consequence of the invariance 
of the Hamiltonian function with respect to a translational shift of the origin 
of coordinates by which the position of the system as a whole is measured. 
Conservation of angular momentum results from the invariance of the 
Hamiltonian function of an isolated system with respect to the angular 
orientation of the coordinate system used to describe it. 

Conservation of spin statistics asserts that the “spin character of a 
closed system cannot change. That is, if a system consists of a number of 
particles of various spins, some integral and some half-integral, the total 
number of particles of half-integral spin must remain either even or odd. 
This conservation law is in some respects a consequence of the conservation 
of angular momentum, since orbital angular momentum can have only 
integral quantum numbers and there fore cannot “carry away’ all of the 
intrinsic angular momentum of a particle. 
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Angular Distribution of Particles 


During the process of nuclear fusion, a deuteron and triton combine 
or fuse to become a compound or excited nucleus. At some point shortly 
(10* 15 to 10' 12 seconds) after the fusion of the two nuclei, the nucleus 
separates into two distinct fragments, a neutron and an alpha particle with a 
total energy of 17.6 MeV. 

These two particles separate along the nuclear symmetry axis 
(References 14,15). The quantum number K, which is the projection of the 
total vector angular momentum quantum number J on the nuclear symmetry 
axis, is a descriptive quantum number describing the passage of a nucleus 
from its excited transition state to the configuration of separated particles. 
The directional dependence of the separation particles resulting from a 
transition state with quantum numbers J, K, and spin quantum number M is 
uniquely determined. 

The quantum numbers J and M (the spin projection of J on a space- 
fixed axis which is usually taken as the direction of motion of the incident 
particle in the center-of-mass coordinate system) are conserved in the entire 
particle separation process. Whereas J and M are fixed throughout the 
various extended shapes on the path to particle separation, no such 
restriction holds for the parameter K . In going from the original compound 
nucleus (with a shape which is the same as the of the ground state) to the 
saddle point or transition state, the nucleus suffers many vibration and 
changes in shape, and redistributes its energy and angular momentum in 
many ways. The K value (or values) of the transition nucleus are, therefore, 
unrelated to the initial K values of the compound nucleus. Once the nucleus 
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reaches the transition state deformation, K is a good quantum number 
beyond this point of the particle separation process. 

The angular momentum coupling scheme for a deformed nucleus and 
the relationship between J, M, and K is schematically illustrated in Figure 2. 



Figure 2. Angular Momentum Coupling Scheme 

The vector J defines the total angular momentum. The quantity M is 
the spin component of the total angular momentum on the space-fixed Z 
axis, which is the direction of the incident particle velocity. The quantity K 
is the component of the total angular momentum along the nuclear 
symmetry axis. The collective rotational angular momentum R is 
perpendicular to the nuclear symmetry axis; thus, K is entirely a property of 
the intrinsic motion. The angle theta represents the angle between the 
nuclear symmetry axis an the space-fixed Z spin. 

The probability of particles being emitted from a transition state with 
quantum numbers J, M, and K at an angle theta is given by: 
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J 

where O represents the probability of emitting particles at angle 

theta into the conical volume defined by the angular increment c* & . The 
normalization is such that the probability integrates to unity for limits 0 and 
pi. The area of the annular ring on a sphere of radius R through which the 
particles are passing is given by the width of the strip R * m fr times the 
circumference of the ring 2flR sin which is 2?R 2 sin i} This annular 
ring area must be divided by the total area of the sphere 4 7R 2 in order to 
give the probability as defined by the above equation. 

j 

The foregoing probability distribution depends on the *Cpdyfu nction 
and is universal in the sense that it is independent of the polar angle, the 
angle of rotation about the symmetry axis, and the moments of inertia. 
Hence the probability distribution depends only on the angle theta between 
the space-fixed and body-fixed axes. The (&) functions are defined by 
the following relation // 


vX r- • r 

5 c-0 L ces 

Z. / \ I 


&)] 




x (y-tf-x)l (j( x *! 


where the sum is over X = 0,1,2,3,... and contains all terms in which no 
negative value appears in the denominator of the sum for any one of the 
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quantities in parentheses. There are a number of symmetry relationships 
among the d’s, for example, 

(->y n (*) 

The angular distribution /f (^' is obtained by dividing the 

probability for emitting particles at angle theta as defined in the above 
equation by sin , 

i/v 7 (*) = fic^+OA] I ^4* 

/M t\ *— 

The normalization of this equation is such that 

f* W^ n (*) = 1 

The angular distribution function C^ 1S directly related to the 

differential particle cross section for a particular channel (J, at angle 

theta by the relation 

J 

ft f j ■» K £) ~ ^‘■y tt (&) 

J/L ** 

where the differential particle cross section is to be expressed in the same 
units per steradian as the total fusion cross section. The factor 2 appears 
instead of the usual 4 It due to the fact that two particles separated by 180 
degrees arise from each fusion event. 

The angular distribution of particles depends on two basic quantities: 
the angular momentum brought in the incident nucleus and the fraction of 
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this angular momentum which is converted into orbital angular momentum 
between the neutron and alpha particle. This fraction is characterized 
through the parameter K, where K is defined as the projection of the angular 
momentum J on the nuclear symmetry axis. These are sufficient to describe 
the emitted particle angular distribution. This is based on the particles 
separating along the nuclear symmetry axis and that K is a good quantum 
number beyond the transition point. The distribution of the angular 
momenta resulting from the fusion system may be controlled by the proper 
choice of spin alignment. 

The density of levels in the transition nucleus with spin J and 
projection of J on the nuclear symmetry axis equal to K is given by the 
approximate relation 

A (^3 / /T) eXfiC ^ y) / 


rot ls the energy tied up in rotation for 

transition state (J,K), and t is the thermodynamic temperature. The 

r r J ' n 

thermodynamic energy available to the nucleus is the quantity £ “ • 

The rotational energy of a nucleus in its saddle point deformation is 

r 7 '"-_ (£/iI x )(?-H‘) f 


where 1^ and I are nuclear moments of inertia about axes perpendicular 
and parallel to the symmetry axis, respectively. Substitution into the above 
equation yields 
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For fixed values of E and t the number of transition levels p (J,K) 

, JL _ »- ' 

4) J 


depends on two quantities, 


ll t ? 


and 


•/*«)- C )] 


If in addition, J is fixed, then the distribution in K becomes 

exp ft /a )/^i C ) * C ' 



This equation is equivalent to a Gaussian K distribution 


sc*)*- e xp {-fry^:) / 


- £ 


/r>x 


where 


y -- -i 


[C ’/I.) - C ‘A> )J 


If the quantity [c 'A.)- C^S) is replaced by 1/I efr , then 

/r/ - 


For a Gaussian K distribution an exact expression for the particle angular 
distribution may be derived by proper weightings of J,M, and K. If the 
target and projectile spin are included, the exact equation for the angular 
distribution utilizing the proper statistical weights is 
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The quantities I D , s, and S are the target spin, projectile spin, and the 
channel spin, respectively. The channel spin S is defined by the relation S = 
I 0 + s . The total angular momentum J is given by the sum of the channel 
spin and the orbital angular momentum, J = S + 1 . The projection of I 0 on 
the space-fixed axis is given by whereas the projection of S ( and J ) on 
the space-fixed axis is M. The quantity in the bracket gives the weighting 
factor for a particular (J,M) combination. This value multiplies the angular 
dependent term for the allowable k states ( K distribution is weighted also) 


of a particular J. This product is summed first over all M values for a 

y 

particular J and finally over all J values. The <£U ff ^fiinction is defined in 

S t Jt t J 

the saddle point deformation equation and the quantities C ^ 0 _ and 

^ Xo, 3 , 3 ' ' 

are Clebsch-Gordan coefficients. 
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